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Abstract 

The effect of edge on wetting and layering transitions of a three-dimensional spin- 
1/2 Ising model is investigated, in the presence of longitudinal and surface magnetic fields, 
using mean field (MF) theory and Monte Carlo (MC) simulations. For T = 0, the ground 
state phase diagram shows that there exist only three allowed transitions, namely: surface 
and bulk transition, surface transition and bulk transition. However, there exist a surface 
intra-layering temperature r£, above which the surface and the intra-layering surface 
transitions occur. While the bulk layering and intra-layering transitions appear above 
an other finite temperature T^{> T£). These surface and bulk intra-layering transitions 
are not seen in the perfect surfaces case. Numerical values of T[ and T^, computed by 
Monte Carlo method are found to be smaller than those obtained using mean field theory. 
However, the results predicted by the two methods become similar, and are exactly those 
given by the ground state phase diagram, for very low temperatures. On the other hand, 
the behavior of the local magnetizations as a function of the external magnetic field, shows 
that the transitions are of the first order type. 

T[ and decrease when increasing the system size and/or the surface magnetic field. In 
particular, reaches the wetting temperature for sufficiently large system sizes. 
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1 Introduction 



Recently, much attention has been directed to study the layering transitions of magnetic solid 
films. Several models have been proposed to study such physical systems. A variety of possible 
phase transitions has been reviewed by de Oliveira and Griffiths [1], Pandit and Wortis [2], 
Pandit et al. [3] and Ebner et al. [4]. Such transitions have been observed in a variety of 
systems including for example ^He [5] and ethylene adsorbed on graphite [6]. Hanke et al. [7] 
show that symmetry breaking fields give rise to nontrivial and long-ranged order parameter 
profiles for critical systems such as fluids, alloys, or magnets confined to wedges. Benyoussef 
and Ez-Zahraouy [8] have studied the layering transitions of a spin-1/2 Ising film using transfer 
matrix and the real space renormalization group method [9]. By using Monte Carlo simula- 
tions on thin Ising films with competing walls. Binder et al. [10], found that occurring phase 
transitions belong to the universality class of the two-dimensional Ising model and found that 
the transition is shifted to a temperature just below the wetting transition of a semi-infinite 
fluid. 

On the other hand, a variation of phase diagrams with the strength of the substrate potential 
in a lattice gas model for multilayer adsorption is studied by Patrykicjcw et al. [11] using 
Monte Carlo simulations and molecular field approximation. Therefore, the effect of finite size 
on such transitions has been studied, in thin film confined between parallel plates or walls, by 
Nakanishi and Fisher [12] using mean field theory and by Bruno et al. [13] taking into account 
the capillary condensation effect. In the framework of the mean field theory, we found in a 
previous work [14], the wetting and layering transitions of three-dimensional spin— 1/2 Ising 
transverse model in the presence of both an external and surface magnetic fields. We showed, in 
our previous work [15], the existence of the layering transitions for a film, with finite thickness 
and infinite surfaces. We have found that the wetting temperature T^, depends weakly on the 
surface corrugation degree for fixed values of the surface magnetic field. 

In addition, by applying a suitable effective interface model at liquid- vapor coexistence, Rejmer 
et al. [16] found a filling transition at which the height of the meniscus becomes macroscop- 
ically large while the planar walls of the wedge far away from its remain non-wet up to the 
wetting transition occurring at T^,. They also showed that the discontinuous filling transition is 
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accompanied by a pre-filling line extending into the vapor phase of the bulk phase diagram and 
describing a transition from a small to a large, but finite, meniscus height. Long-ranged order 
parameter profiles, for critical systems such as fluids, alloys, or magnets confined to wedges, 
have been studied by Hanke et al. [7]. They also discuss the properties of the corresponding 
universal scaling functions of the order parameter profile and the two-point correlation function, 
and determine the critical exponents for the normal transition. 

The aim of this work is to study the effect of edge surfaces on the wetting and layering 
transitions of a three-dimensional spin— 1/2 Ising model in the presence of external and sur- 
face magnetic fields, using mean field theory (MF) and Monte Carlo (MC) simulations. It is 
found that there exist a surface intra-layering temperature Tf, above which the surface and the 
intra-layering surface transitions occur. While the bulk layering and intra-layering transitions 
appear above an other finite temperature T|(> T£). Such results were not seen neither in the 
perfect surfaces case [2,3,8,9,11,14,15] nor in the continuous model case [7,16]. The outhne of 
this paper is as follows. In the following section we describe the model and the formulations 
used. Section 3 is devoted to results and discussions. 

2 Model and method 

The system we are studying. Fig. 1, is formed with N layers. Each layer '/c' [k — 1,2, N), 
is formed with two perpendicular perfect planes and contains 2(A^ — k) + 1 spin chains which 
are infinite in the y— direction. A uniform surface magnetic field Hg is applied on the planes 
{x, y,z = 1) , {x = 1, y, z) of the layer k = 1. However, the system is invariant by translation, 
in the y— direction, so the coordinate y will be dropped in the following. An external field H 
is applied on the global system. The Hamiltonian governing the system, can be written as 

= - X! Jij^i^j ~ ^HiSi (1) 

<i,j> i 

where. Si — ±1 are spin— 1/2 Ising random variables, and Jij — J are the exchange interactions 
assumed to be constant. The system is assumed to be infinite in the direction y, so the variable 
y will be cancelled in all the following. is the total magnetic field applied on a spin located 
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on a site 'i' of coordinates {x, z). defined by: 



Hi = H{x, z) = < 



H + Hs for X = = 1,...,A^ 

H + Hs for z = l,x = l,...,N (2) 
H elsewhere 



H and Hg are, respectively, the external and surface magnetic fields both applied in the 
direction. 

Using the mean field theory, we compute the magnetizations and the free energy of such a 
system. However, the magnetizations can be expressed as: 

m(x, z) — tanh(/3(2m(x, z)-\-m{x, z-\-l)-\-m{x, z — l)+m(x + l, z)-\-m{x — l, z)+H{x, z))) (3) 

where /3 = 1/(A;bT), ks denotes the Boltzmann constant and T the absolute temperature. 
Eq.(3) of the magnetizations is the Callen identity with an external field H{x,z) so that the 
total free energy of the system can be derived as : 

F[m{x, z)] = Y.x,z{'^[ ^~"'t'''^ ^"9{^ - z)) + i±^^^^Lo5((l + m{x, z))] 

— ^m{x, z) [2m(a;, z) + m{x, z + 1) + m{x, z — 1) + m{x + 1,2;)+ m{x — l,z)] (4) 

— m(x, z)H{x, z)}. 

The system is with free boundary conditions. 



3 Results and discussion 

The notation (Fl^) with p = 0, 1, 2, TV and g = 0, 1, 3, 5, 2{N - p) - 1, will be used to 
denote that the first p layers and the q first spin chain of the layer p + lare in a magnetic 
state "up"; while the remaining N — (p + l) layers are in the state "down". In particular, the 
notation 1^ (resp. O^) will be used to denote a configuration with positive magnetization for 
all layer spins (resp. negative magnetization for all layer spins) of the system. Fig. 2 illustrates 
an example of these notations for N — A. 

The ground state energy of a given configuration is calculated exactly from the Hamiltonian 
(1). The transition from a configuration (l^l^) to an other configuration (l^'l^/) is obtained 
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by the equality of their energies. It is found that the ground state transitions, illustrated by 
Fig. 3, occur according to linear forms H/ J — a -\- b{Hs/ J). We found three ground state 
transitions: 

• The surface transition ^ (IHq) located at H/J = -HjJ + 2{N - 1)/{2N - 1). 

• The bulk transition (l^lo) ^ 1^ occurring at H/J = 2/(1 — N). 

• The bulk and surface transition ^ 1^ locahsed at H/J = ((1 - 2N)/{N^))Hs/ J. It 
is clear that for large values of the system size N, (N — > oo), the bulk and surface, and 
the bulk transitions are shifted to H/J = 0. Whereas the surface transition is located at 
H/J = —Hs/ J + 1. This is in good agreement with the results we have established in our 
previous works [13,15], for the infinite perfect surfaces case. 

Although the established equation are valid for an arbitrary system size N , numerical results 
are given for two system sizes iV = 4 (thin film) and = 20 (thick film) spins in both directions 
X and z. In order to examine the effect of temperature on wetting and layering transitions, we 
plot in Fig. 4 the corresponding phase diagram, for a fixed value of the surface magnetic field 
Hg, by using the mean field (MF) theory. Fig. 4 is plotted for Hs/ J = 1.0 and a small system 
size N — 4. From the ground state phase diagram (see Fig. 3), the situation corresponding 
to Fig. 4 is located in the region: bulk transition. Indeed, for very small temperature values, 
the only transition seen is 1^. While the increasing temperature leads to the intra- 

layering surface transitions and the intra-layering bulk transitions as well as the inter-layering 
transitions. In the case of = 4, Fig. 2 shows an example of the intra-layering transitions, 
which correspond to the change of the spin chains from "down" to "up" state inside the layer 
(in surface and bulk). One can note that these transitions arc due to the geometry of the 
surfaces, independently on the system size A". Indeed, the intra-layering surface transitions 
are: O' ^ {I'h), (iH,) ^ (iH,), (iH,) ^ (1^,), (iH,) ^ (l^lo) and (I^Iq) ^ l'; 
while the intra-layering bulk transitions are ^ (l^li); (l^li) ^ (l^ls); (l^ls) ^ (l^lo), 
(l^lo) ^ (l^li), (l^li) ^ (l^lo) and (iHo) ^ 1" (sec Fig. 2). 

The profile of the magnetizations m(l, 1), m(l,2), m(2,2), m(2, 3) are given in Figs. 5a and 
5b. From these figures, it is clear that the configuration (l^^li) is seen ( Fig. 5a) before the 
configuration (l^ls). The same argument is valid for the configurations (l^li) and (l^ls): Fig. 
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5b. 

On the other hand, the Monte Carlo calculations are performed on a system with a total number 
of spins 4 X 4 X 100 = 1600 (N — A and % = 100 spins in the y-direction) . The phase diagram 
obtained using the Monte Carlo (MC) method is illustrated by Fig. 6 for the same surface 
magnetic field value used in the mean field method Hg/ J = 1.0 (Fig. 4). Comparing these 
figures. (Fig. 4 and Fig. 6), plotted for the same surface magnetic field value, it is seen that the 
layering transitions are found by the (MC) method, when increasing the temperature, before 
these same transitions can be seen by the (MF) method. But the global transitions are found 
by the two methods, and qualitatively the phase diagrams obtained exhibit similar topologies. 
It is worth to note that for very low temperature, (MF) and the (MC) are in good agreement 
with the ground state phase diagram. 

In order to outline the above results, we plot the local magnetizations m(l, 1), m(l, 2), m(2, 2) 
and m(2, 3) in Figs. 7a and 7b, as a function of the reduced bulk field H/ J. When increasing 
the external field H/ J the corners transits (Figs. 7a and 7b) before the other intra-layering 
transitions. Such behavior was not obtained in the perfect surface case model. However, the 
other layer-by-layer transitions occur when increasing the external field at specific temperature 
values. As one can expects, the external field values needed to make arising the intra-layering 
and inter-layering transitions increase with increasing the order of the layer counted from the 
surface k = 1. This is qualitatively in a good agreement with our previous works [13-15]. 
The phase diagram of the intra-layering and layer-by-layer transitions, in the plane {H/ J, T/ J) 
for a larger system size with = 20 and Hg/ J — 1.0, is plotted in Fig. 8, by using the mean 
field method. Comparing this figure with Fig. 4, it is seen that the increasing system size efi'ect 
is to decrease the intra and inter layering temperature values. The bulk layering transitions are 
shifted to higher external magnetic field values, for a fixed surface magnetic field value. 
To complete this study, we have investigated in Fig. 9 the surface, T|, and bulk intra-layering 
temperature profiles as a function of the system size N for two surface magnetic field values 
Hs/J = 1.0 and Hg/J = 0.9. It is found that (Fig. 9a) (as well as T|, Fig. 9b) decreases 
when increasing the system size. For a fixed system size N , these intra-layering temperatures 
decrease when increasing the surface magnetic field values. Moreover, reaches the wetting 
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temperature for sufficiently large values of the system size. 

To conclude this study, we note that the presence of an edge (angle 7r/2) decreases the wetting 
temperature compared with the perfect case (angle tt), see ref. [14]. This finding is in a good 
agreement the results obtained by Rejmer et al. [16]. 

4 Conclusion 

The intra and inter-layering transitions of a three-dimensional spins Ising model with edge 
surfaces are studied within the mean field theory and Monte Carlo simulations. At T = 0, 
the only layering transitions, occurring under the effect of an external and a surface magnetic 
fields, are: surface transition, bulk transition and surface to bulk transition. When increasing 
the temperature, a succession of intra-laycring transitions, absent in the case of perfect surfaces 
[13,15] and the continuous model [7,16], are found. The layering temperatures obtained by the 
mean field theory are higher than those predicted by MC method. The two methods lead to 
similar results for very low temperatures. The effect of the edge on the behavior of the local 
magnetizations m{x, z), as a function of the reduced bulk field H/ J, is also investigated. We 
have analyzed the effect the system size on the surface and bulk layering temperatures as well 
as on the wetting temperature. 
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Figure Captions 



Figure 1: Geometry of the system formed with two surfaces {x, y,z — 1) and {x — l,y, z) with 
N spins in both the x and 2;— directions. The system is infinite in the direction. A uniform 
surface magnetic field Hs is apphed on the planes {z = l,x,y) and {x = l,y,z). An external 
magnetic field H is apphed to the global system. 

Figure 2: Sketch of different possible configurations for a system with N — A. Symbols (o) 
and (•) correspond to spin "down" and spin "up" respectively. The notation (Pig) where 
p = 0, 1, 2, N and q = 0,1, 3, 5, 2{N — p) — 1, will be used to denote that the first p layers 
and the q first spins of the layer p + 1 are with positive magnetizations, while the remaining 
N — [p+l) layers are with negative magnetizations. In particular, the notation 1-^ (resp. O^) 
will denote a configuration in the state "up" for all layer spins (resp. state "down" for all layer 
spins) of the system. 

Figure 3: The ground state phase diagram in the {H/ J, Hs/ J) plane. There exist only three 
transitions, namely: the surface layering transition ^ (l^lo) , the surface to bulk layering 
transition (l^lo) 1^ and the bulk layering transition <-> 1^. 

Figure 4: Layer- by-layer and intra-layering transitions, in the plane [H / J,T / J), using the 
mean field method for a system with N — A and Hgj J — 1.0. 

Figure 5: Magnetization profiles as a function of the reduced bulk magnetic field HjJ for 
= 4 and EsjJ = 1.0, of m(l, 1) and m(l,2) at T/J = 3.95 (a), m(2,2) and m(2,3) at 
TjJ = 4.75 (b), by the mean field method. 

Figure 6: Phase diagram of the intra-layering and layer-by-layer transitions, in the plane 
{H/ J,T / J), using the Monte Carlo simulations for a system with N — A,ny — 100 and 
HJJ = 1.0. 
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Figure 7: Magnetization profiles as a function of the reduced bulk magnetic field H/ J for 
= 4 and Hs/J = 1.0, of m(l, 1) and m(l,2) at T/J = 0.75 (a); m(2, 2) and m(2,3) at 
T/J =1.0 (b), using the Monte Carlo method. 

Figure 8: Phase diagram of the intra-layering and layer-by-layer transitions, in the plane 
{H/ J, T/J), using the mean field method for a system with N — and Hg/ J — 1.0. 

Figure 9: Surface (a) and bulk (b) layering temperature profiles as a function of the system 
size N for two surface magnetic field values Hs/J =1.0 and Hs/J = 0.9, by using the mean 
field method. 
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